Charging of a quantum dot coupled to Luttinger liquid leads 
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Luttinger liquid behavior of one-dimensional correlated electron systems is characterized by power- 
law scaling of a variety of physical observables with exponents determined by a single interaction 
dependent parameter K. We suggest a setup to study Luttinger liquid behavior in quantum wires 
which allows to determine K from two independent measurements: transport through a quantum 
dot embedded in the wire and the charge on the dot. Consistency of the two K's for a single probe 
would provide strong experimental evidence for the Luttinger liquid paradigm. 
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I. INTRODUCTION 

Theoretically it is well established that the two-particle 
interaction U in metallic, one-dimensional (Id) electron 
systems leads to Luttinger liquid (LL) physics.— One of 
the characterizing properties is the power-law scaling of 
a variety of physical observables as functions of exter- 
nal parameters (e.g. the temperature T) with exponents 
which can be expressed in terms of a single, interaction 
dependent LL parameter K < 1 (for repulsive interac- 
tions, with K = 1 for U = 0). However, there are only a 
few experiments on quasi Id systems which reveal clear 
indications of LL behavior. Even in these rare examples 
mostly only a single observable as a function of a single 
control parameter was measured. In this case it is very 
difficult to convincingly exclude any other source for the 
observed power- law scaling than LL physics. The situa- 
tion becomes more complex as even for the same type of 
quantum wire (e.g. semi-conductor heterostructures, car- 
bon nanotube, or a row of atoms on a surface) K might 
vary from probe to probe since it depends not only on U 
but also on other details such as the band structure and 
filling^ 

A more direct evidence for LL behavior could be 
achieved in the following way. Using a single probe one 
should measure two observables for which LL theory pre- 
dicts power-law scaling with different exponents (3\{K) 
and 02 (K) - in the optimal case even as functions of two 
different control parameters. If the two exponents turn 
out to be consistent, that is K{(3\) « K^fa), strong evi- 
dence for LL physics is achieved. A step in this direction 
is the linear conductance G(T) measurement by Yao et 
al£ across an impurity free part of a metallic single-wall 
carbon nanotube as well as across a part of the same tube 
containing a single kink (impurity). 

We propose the setup sketched in Fig. [TJ in which one 
can measure G through a quantum dot (QD) embedded 
in a Id wire as a function of T and, for the same probe, 
the charge n accumulated on the dot as a function of the 
dot level position varied by an external gate voltage V g . 
The QD is formed by two high barriers within the Id 
system, e.g. realized by additional gates. The charge is 
detected by the current running through a nearby quan- 
tum point contact (QPC). While transport of Id corre- 
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FIG. 1: (Color online) Sketch of the setup to measure trans- 
port through a quantum dot as well as the charge on the dot. 



lated electrons through double barriers has extensively 
been studied theoretically in recent years,— less is known 
about the charging of a small QD coupled to two LL 
leads. 4 Here we investigate in detail how n(V g ) is affected 
by LL physics using two approaches. First we consider a 
field-theoretical, effective low-energy model, the infinite 
Luttinger model, 1 and perturbation theory in the dot-LL 
coupling r. This can be done for arbitrary < K < 1, 
but is restricted to small T, In a complementary, second 
step we study an interacting microscopic lattice model 
of finite length coupled to noninteracting leads, a model 
being closer to experimental setups. To treat the cor- 
relations we use the functional renormalization group 
(fRG)£ This method can be applied for arbitrary T, but 
is restricted to small 1 — K. Both approaches lead to 
consistent results and we show that n(V g ) is governed 
by power-law scaling, which should be detectable in the 
suggested setup. We consider a dot with a large level 
spacing such that only a single level matters. Further- 
more, we mainly consider spinless fermions and suppress 
the Kondo effect. Experimentally this can be achieved 
by a magnetic field lifting the spin degeneracy of the dot 
level or by measuring at T > Tk, with Tk being the 
Kondo temperature. For transport through a dot show- 
ing the Kondo effect coupled to LL leads, see Ref. @- 

From the linear conductance G(V g ,T) through a dou- 
ble barrier K can be extracted in several ways, some of 
them restricted to certain regimes of K values or sym- 
metric barriers.— To be as general as possible we here 
present a prediction which holds for all < K < 1 , sym- 
metric as well as asymmetric barriers, and which does not 
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require any other fine tuning of parameters. For a fixed 
gate voltage away from resonance, which we assume to be 
at V g = 0, one finds G ~ max{T, ^} 2 ( 1 / K_1 ) at asymp- 
totically small scales.-' 7 Here 8 denotes an energy scale 
~ l/N, with N being the length of the LL wire, which is 
eventually coupled to noninteracting leads. 

In an important work Furusaki and Matveev analyzed 
n(V g ) for strongly interacting systems with K < 1/2 
within the infinite Luttinger model using perturbation 
theory in T and the mapping to related problems^ They 
showed that for sufficiently small T, n(V g ) is discontinu- 
ous at V g = 0. For 1/3 < K < 1/2 the finite V g behavior 
adjacent to the jump shows scaling with the exponent 
1/K — 2, while for even smaller K the deviations are lin- 
ear in V g . The perturbation theory in T for the Green 
function - not for the self-energy, as used by us - breaks 
down for 1/2 < K < 1. In an attempt to investigate LLs 
characterized by such if's a numerical method was used 
for systems of up to 150 sites in Ref. |9|. The authors con- 
cluded that n(V g ) is continuous and does not show LL 
physics. Below we confirm the first statement but show 
that the second is incorrect as finite size corrections com- 
pletely mask the power-law behavior. 



The dot charge is 



n(V g ) 



du: p d {uj) , 
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with the chemical potential fi = 0. Because of the 
particle-hole symmetry it obeys n(V g ) = 1 — n{— V g ) and 
from now on we focus on 7 9 > 0. In contrast to the 
perturbation theory in T for the dot Green function it- 
self used in Ref. |8| which is restricted to K < 1/2, our 
approach can be applied for all < K < 1. 

Based on Eqs. ©-(Ill) the leading small V g behav- 
ior of 1/2 — n(V g ) can be determined analytically. For 
1/2 < K < 1, n(V g ) is a continuous function with 
n(V g = 0) = 1/2. This implies that the width w over 
which n(V g ) changes from 1 to is finite.— The func- 
tion n(V g ) contains regular terms proportional to V g 2l+1 , 
with I G No, as well as anomalous terms with expo- 
nents containing K. The leading anomalous term is 



V, 



(2K-1)/(1-K) 



Depending on K either the linear term 



or the anomalous term dominates. A special situation is 
reached at K = 2/3, where logarithmic corrections ap- 
pear. The leading V g dependence is given by 



II. PERTURBATION THEORY IN THE 
LEVEL-LEAD COUPLING FOR THE 
SEMI-INFINITE LUTTINGER MODEL 

We first consider a QD coupled to two LLs via tunnel 
barriers with hopping amplitudes ti/ r . For simplicity the 
LLs are assumed to be equal and described by the semi- 
infinite Luttinger model 1 (with an open boundary on the 
side coupled to the dot). To leading order in T = tf + 1 2 
the dot self-energy is given by T,d(z) = TQ(z), with the 
single-particle Green function Q of the disconnected semi- 
infinite LL at the boundary. The low-energy behavior of 
the imaginary part oiQ for z = u+iO, that is the spectral 
function p, is known exactly from bosonizationiiS It is 
given by p(uj) ~ |oj| 1//a_1 . To be specific we assume that 
p(u>) has support [— u> c ,uj c ] 
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It is then straight forward to compute ReQ R (u>) by 
Hilbert transformation. The leading behavior at 
\u)/lo c \ <C 1 is given by 
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Using the Dyson equation the dot spectral function pd 
follows from the perturbative as 



Pd{uj) = 



r P (w) 



[w — V g — TRe£«H] 2 + brrp(w)] 2 



(3) 
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for | < K < 1 
for K = | 



(2K-l)/(l-K) 
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for i < K < 



At K = 1/2, n(Vg) is still continuous and for V g \ 
approaches 1/2 with corrections ~ 1/| ln(V g /uj c )\. 

For K < 1/2 and small T, n(V g ) shows a jump at 
V g = 0, that is ]im Vg \o n(V g ) = A < 1/2. In this 
regime our perturbative approach for the self-energy, 
which guarantees the correct analytical structure of the 
dot Green function, gives the same results as the pertur- 
bation theory for the Green function itself used in Ref. 8. 
This follows from two observations. According to Eq. @ 
the real part of Q R becomes linear at small u> and can 
thus be absorbed in the first term in the denominator of 
Eq. ([3]). In addition, for small V g the contribution of p 
in the denominator of Eq. ([3]) can be neglected compared 
to the term linear in u>. For the small V g analysis and 
to leading order in T the integrand in Eq. ^ becomes 
equivalent to the one obtained in Ref. [H 
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The jump at V g \ is given by A = T/[(2 - ^K)lj 2 c ] 
which is nonuniversal as it depends on the cutoff u> c . Ev- 
idently, for K close to 1/2 this expression only holds for 
sufficiently small T. In Ref. 8 it is argued that increas- 
ing r beyond the perturbative regi me A decreases, ap- 
proaches the minimal value Ao = y K/2 at a certain T , 
and for T > Tq, n becomes a continuous function of V g 
even for K < 1/2. The finite V g corrections of n for small 



3 



r are given by 



A - n(V g ) ~ { Yr. 



for | < K < \ 




(J) 



forO < K < ± 



These results show that for 1/3 < K < 2/3, that is for 
sufficiently strong, but not too strong interactions, the 
LL parameter K can be extracted from a measurement 
of n(V g ) for gate voltages close to the resonance value. 

A second way to extract the LL parameter in the 
regime in which n(V g ) is continuous, that is for 1/2 < 
K < 1, is given by the T dependence of the charac- 
teristic width w over which the charge changes from 1 
to 0. In particular, this includes weak interactions with 
2/3 < K < 1 for which 1/2 - n(V g ) itself is linear in V g 
and cannot directly be used to determine K . The width 
can e.g. be defined by w = 2V g with n(V g ) = 1/4. In ex- 
perimental setups in which the two barriers are realized 
by gates, V can be tuned by varying the applied voltages 
and w(T) can be extracted. For r — > 0, w(T) follows 
from Eq. ^ and scales as 



w(T) 



r 



K/(2K-l) 



for 1/2 < K< 1 



(8) 



On first glance the appearance of an anomalous exponent 
in w might be at odds with the linear V g dependence of 
1/2 - n(Vg) for 2/3 < K < 1. In fact, both results are 
consistent as the regime over which n(V g ) goes linearly 
through 1/2 around V g ~ shrinks with decreasing F and 
decreasing K. To experimentally observe the predicted 
power-law scaling the temperature has to be sufficiently 
smaller than the width w. 

In the absence of the Kondo effect (see above) including 
the spin degree of freedom does not lead to new physics. 
The perturbative analysis can be repeated after replacing 
the exponent 1 /K — 1 in the spectral function of Eq. (TT|) 
by the exponent for LLs with spin (1/K — l)/2>i 



III. WEAK TO INTERMEDIATE 
INTERACTIONS IN A MICROSCOPIC LATTICE 
MODEL 

We next replace the LLs described by the semi-infinite 
Luttinger model by the microscopic lattice model with 
nearest-neighbor hopping t > and nearest-neighbor in- 
teraction U . On both sides of the QD the LLs are as- 
sumed to be finite, each having sa N/2 sites, and adiabat- 
ically coupled to noninteracting Id tight-binding leads£ 
The interaction is treated by an approximation scheme 
that is based on the fRG and which was shown to be re- 
liable for weak to intermediate interactions^ In contrast 
to the perturbation theory in T, which is restricted to 
small r, this method can be applied for all T and is thus 
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FIG. 2: (Color online) Main part: Functional RG results for 
the dot spectral function as a function of energy for U/t = 0.5, 
ti/t = t r /t= VoJ, V g /t = l, and length of the LL wire N = 
10 2 (dashed dotted line), N = 10 4 (dashed line), N = 10 5 
(solid line). Note the log-scale of the j/-axis. Lower inset: 
Zoom-in of pa around u> = (linear-scale of y-axis). Upper 
inset: Dot occupancy as a function of the gate voltage for the 
same parameters as in the main part and -/V = 10 5 . 



complementary to the above approach. The Hamiltonian 
is given by 



H = —t 



j = -oo 
JV-1 



— — OO 

(*' - *)4 c «-i - (*r - t)c] d+1 c 3d - H.c. (9) 



c]cj. 



in standard second quantized notation, with nj 
To prevent any backscattering from the contacts to the 
noninteracting leads around j ~ 1 and N the interaction 
is turned on and off smoothly over a few lattice sites, with 
a bulk value U, as described in Ref. 5. The dot is located 
at lattice site jd somewhere close to N/2 (the results are 
insensitive to the exact position). The prime at the sum 
in the second line indicates, that the interaction across 
the barriers defining the QD is set to zero. We also stud- 
ied the case in which the interaction on these bonds takes 
the bulk value U and found that our conclusions are valid 
also for this setup. We focus on half-filling of the band. 
In this case the bulk model is a LL for \U\ < 2t and a 
closed expression for K in terms of the model parameters 
can be giver* 1 ^ K^ 1 = 2arccos [—U/(2t)] /it. 

Within the fRG one introduces an energy cutoff A into 
the noninteracting propagator. Taking the derivative of 
the generating functional of the one-particle irreducible 
vertices with respect to A and neglecting higher order 
corrections one derives a set of O(N) coupled differential 
equations for the A-flow of the self-energy and a renor- 
malized nearest-neighbor interaction. It can be solved 
numerically for up to 10 7 sites, resulting in an approx- 
imate expression for the dot Green function. This ap- 
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FIG. 3: (Color online) Functional RG results for the width 
over which n(V g ) changes from 1 to as a function of the 
dot-LL coupling T for TV = 10 5 and different U. 
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n(V g ) by a Fermi liquid form. 

The LL suppression of pd around tu = will manifest 
itself also in the charging of the dot. To illustrate this 
we confirm the prediction of Eq. ((8]) for the T depen- 
dence of w. We extract vjH from n(V g ) (for an example 
of n{V g ) see the inset of Fig. [5]) for N — 10 5 , a variety 
of r (for simplicity assuming symmetric barriers), and 
different U. The results for w(T) are shown in Fig. [3] 
on a log-log scale. At small T, w shows power-law scal- 
ing. In Fig. 0] the exponent as a function of U, obtained 
by fitting the data of Fig. [3] (and additional data sets), 
is compared to K/(2K — 1) determined in perturbation 
theory in T [see Eq. (JSj)] . We used the exact relation 
between K and U mentioned above. The results agree 
quite well for < U/t < 1/2. For larger U higher or- 
der corrections neglected in our truncated fRG scheme 
become important. For sufficiently large U the exponent 
K/(2K— 1) becomes large (it diverges for K \ 1/2) and 
should experimentally be clearly distinguishable from the 
noninteracting value 1. 
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FIG. 4: (Color online) Circles: Scaling exponent of the width 
w(T) over which n(V g ) changes from 1 to extracted from the 
fRG data for the microscopic model. Solid line: The exponent 
K/(2K-1) [see Eq. ©], with K = K(U), as obtained within 
the Luttinger model. 

proach is described in detail in Ref. [H. 

From the Green function the spectral function pd and 
thus the charge on the dot [see Eq. (|4|)] can be computed. 
In Fig.f2]pd(w) is shown for U = 0.5, symmetric barriers 
tt/t = t r /t = x/OT, Vg/t = 1, and different N (note 
the log-scale of the y-axis in the main part). The upper 
inset shows n(V g ) for N = 10 5 . On the scale of the 
plot n(V g ) does not change if one further increases N. 
The dominating feature of pd is the Lorentzian-like peak 
at u ~ V g . Although a fermion occupying the dot is 
assumed to be noninteracting with the fermions in the 
leads, increasing N the coupling to the LL wires clearly 
leads to a power-law suppression pd{u) ~ uj 1 ^^ 1 close 
to oj = 0, as also given by the perturbative expression 
Eq. ([3]). The lower inset of Fig. [5] shows a zoom- in of the 
dip region. The absence of this LL feature at small N 
of order 100 explains why in Ref. |9| it was possible to fit 



IV. SUMMARY 

Using two different models and methods we have inves- 
tigated the charge n(V g ) accumulated on a QD coupled 
to two LL wires when the dot level position is varied by 
an external gate voltage. Depending on the strength of 
the two-particle interaction U , LL physics manifests it- 
self in power-law scaling of n(V g ) close to the resonance 
at V g = and the width w(T) over which n(V g ) changes 
from 1 to 0. The corresponding exponents can be ex- 
pressed in terms of the LL parameter K . We proposed a 
setup which simultaneously allows to measure n(V g ), and 
thus w(T), as well as the temperature dependence of the 
linear conductance G(T) through the QD. Off-resonance 
the latter is also governed by power-law scaling with an 
exponent which can be expressed in terms of K. Consis- 
tency of the extracted K's would provide strong evidence 
for the experimental observation of LL physics. 
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